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<^ ABSTRACT 

We present results from integral field spectroscopy of a field located near the Trapezium 
Cluster using the Potsdam Multi- Aperture Spectrograph (PMAS). The observed field contains 
a variety of morphological structures: five externally ionized protoplanetary discs (also known 
as proplyds), the high- velocity jet HH 514 and a bowshock. Spatial distribution maps are 
Q obtained for diff'erent emission line fluxes, the c(tlj3) extinction coefficient, electron densities 

^ and temperatures, ionic abundances of different ions from collisionally excited lines (CELs), 

C/^ C^^ and O^^ abundances from recombination lines (RLs) and the abundance discrepancy factor 

of 0^+, ADF(0^+). We distinguish the three most prominent proplyds (177-341, 170-337 
and 170-334) and analyse their impact on the spatial distributions of the above mentioned 
quantities. We find that coUisional de-excitation has a major influence on the line fluxes in 
^ the proplyds. If this is not properly accounted for then physical conditions deduced from 

commonly used line ratios will be in error, leading to unreliable chemical abundances for these 
objects. We obtain the intrinsic emission of the proplyds 177-341, 170-337 and 170-334 by a 
direct subtraction of the background emission, though the last two present some background 
contamination due to their small sizes. A detailed analysis of 177-341 spectra making use 
V/^ of suitable density diagnostics reveals the presence of high-density gas (3.8 x 10^ cm"^) in 

contrast to the typical values observed in the background gas of the nebula (3800 cm"^). We 
also explore how the background subtraction could be aff'ected by the possible opacity of 
^-H the proplyd and its eff'ect on the derivation of physical conditions and chemical abundances 

^ of the proplyd 177-341. We construct a physical model for the proplyd 177-341 finding a 

• i-H good agreement between the predicted and observed line ratios. Finally, we find that the use 

of reliable physical conditions returns an ADF(0^^) about zero for the intrinsic spectra of 
V-j 177-341, while the background emission presents the typical ADF(0^^) observed in the Orion 

Nebula (0.16 ± 0.1 1 dex). We conclude that the presence of high-density ionized gas is severely 
aff'ecting the abundances determined from CELs and, therefore, those from RLs should be 
considered as a better approximation to the true abundances. 

Key words: ISM: abundances; HII regions; - ISM: individual objects: (170-334, 170-337, 177- 
341, Orion Nebula); stars: pre-main- sequence; protostars; planets and satellites: protoplanetary 
discs 
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1 INTRODUCTION 

Circumstellar discs are very common around young stars and it is 
well established that they are an integral part of the star-formation 
and planet-formation processes ( [Williams & Cieza|2011| l. Rotation- 
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ally supported and geometrically thin discs of dense molecular gas 
and dust persist throughout the classical T Tauri stage of low-mass 
pre-main sequence stellar evolution ( Beckwith et al. 11990 ), with 
a duration of several million years. When such discs are found in 
or near an H ii region their detection, structure and evolution are 
strongly influenced by the ultraviolet radiation from hot massive 
stars that dissociates and ionizes the local interstellar medium. The 
currently accepted model ( Sutherland 1997 , Henney & Arthur 1998 , 
Johnstone et a l.4998','Richling & Yorke,1998HStorzer & Hollen-| 



bach|1999[ [Nguyen et al. 2002 , Va sconcelos et al.|2011| l is that far 
ultraviolet photons heat the surface of the disc, forming a warm 
outflowing envelope of neutral gas that becomes photoionized after 
expanding to a few times the outer radius of the disc. A schematic 
illustration of this photoevaporating flow model is shown in Fig. 5 of 
[Henney & O 'Dell ( 1999 ). This special class of protoplanetary disc 
is today known as proplyd, a term that was coined by O'Dell et al.[ 
(T993 ) to describe the tear-drop shaped objects observed in the first 
imaging studies of Orion Nebula using the Hubble Space Telescope 
{HST). Its proximity (436 + 20 pc;'0'Dell & Henney"2008') and high 
surface brightness make the Orion Nebula an important laboratory 
for the detailed study of phenomena associated with star-formation 
and the early stages of stellar evolution. Indeed, the majority of 
currently known proplyds 200 objects, [Ricci et al.|2008| l have 
been discovered in the Orion Nebula, although an increasing number 
are now being detected in other H ii regions ( Stapelfeldt et al.|1997| 



Stecklum et al. 1998 , Brandner et al. 2000, Smith et al. 2003 [Yusef-J 



Zadeh et al.|2005,pe Marco et al. 20 06, Wrigh t et al.,20T2l l. In 



the central part of the Orion Nebula, the Huygens region, the main 
features of the proplyds are the brightest zone of the ionization front, 
known as cusp, and the dusty tails of dense gas facing away from 
6^ Ori C, the principal ionizing star of the Orion Nebula. However, 
what we see depends strongly on the orientation and localization of 
the proplyds, appearing for example as dark silhouettes against the 
bright nebular background when found at a distance from 6^ Ori C 
or within the foreground Orion Veil ( ^McCaughrean & 0'Dell|1996| 
[O'Dell &Wong|1996| ). 

The first ob servations of the Orion proplyds were made by 
Laques & Vidall ([l979 ) who detected six of them around 6^ Ori C 



as unresolved nebular condensations (given LV designations) with 
strong emission in Ha, H/? and [O iii] /15007. From the Hof surface 
brightness of the condensations, [Laques & VidaT| ( | 1979 1 determined 
electron densities ranging from 1.5 X 10^ to 4.1 x 10^ cm~^, al- 
though these were only lower limits since the objects were not fully 
resolved. The correct interpretation of the objects as externally ion- 
ized circumstellar discs was initially proposed on the basis of radio 
interferometric mapping ( Churchwell et al.|1987] l and ground-based 
optical spectroscopy ( Meaburn 1988 ), and this was well confirmed 
by subsequent high-resolution imaging with the HST ( [O'Dell et al.| 
|1993[[O^Dell & Wen|1994|l and ground-based adaptive optics sys- 
tems ( [McCuUough et al.||1995| ). The distribution, evolution and 
structure of proplyds and their embedded discs have been studied 
via imaging at optical, infrared and radio wavelengths (e.g. Church] 



well et al. 


1987|[Chen et al.|1998||Bally et al.,1998.,2000„Lada et^ 


2000 
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bberto et al.|2002| [Smith et al.|2005| Eisner & Carpenter[ 


2006 


[Eisner et al. [[2008V Many spectroscopic studies have also 



been carried out, with the earliest being devoted to the kinematical 
properties of the LV knots through analyzing the [O iii] /15007 emis- 
sion profile with high-resolution echelle spectroscopy (^Meaburn[ 
[T9881 [Meaburn et al.|1993T[Massey & Meaburn] 1 993] [1995|l and 
comparing these results with the first theoretical models (' Henney[ 
[et al. 1 1997) . [Bally et al.[ ( [1998 ) carried out intermediate resolution 
spectroscopy with EOS at HST in the ultraviolet range 1 150-3300 



A for proplyds 15 9-350 and 158-327, detecting Cm], Cii], On] 
and Mg ii emission. [Henney & 0'Dell[ ( p^ 99 ) presented a detailed 
analysis of four proplyds (170-337, 177-341, 182-413, and 244- 
440) performing high-resolution echelle spectroscopy with the 10m 
Keck Telescope. This study analysed the spatial profiles of several 
emission lines. Yip, [Oiii] ^4959, Hei i5676, [Oi] i6300, [Siii] 
/16312, Hof, [Nil] /16583 and [S ii] /16731, as well as the comparison 
with theoretical models providing direct measurements of diff'erent 
proplyd properties such as size, ionized density and outflow velocity. 
Particularly, [Henney & 0'DeIT| ( [1999| l estimated peak densities in 
the cusps of about 10^ cm~^. 

The largest number of detailed spectroscopic studies to date 
have concentrated on the proplyd LV2 (167-317), which is the 
most luminous of the inne r group of proplyds (separation from 
6*^ Ori C is < 10'').|Henney et al. (2002) carried out high-resolution 
(0.05'' X 2 km s~^) NUV long-slit spectroscopy with the STIS spec- 
trograph, determining a peak electron density of ^ 3 x 10^ cm~^ 
in the cusp of this proplyd from the [C iii] .11907/C iii] X1909 line 
ratio. Vasconcel os et al.[p005| l performed ground-based integral 
field spectroscopy of LV2 with CMOS on the 8m Gemini Tele- 
scope, detecting multiple forbidden and permitted emission lines 
in the spectral range from 5500 to 7600 A and analysing the high- 
velocity emission associated with the jet (HH 514; Meaburn 1988) 
that emerges from this proplyd. More recently, Tsamis et al. ( 2 01l| 
and Tsa mis & Walsh ( 201 1 ) have obtained the deepest spectroscopic 
observations of LV2 and its redshifted jet with the integral field spec- 
trograph FLAMES on the 8m VLT in the spectral range from 3000 
to 7200 A. The FLAMES data were complemented with analysis 
of the ultraviolet and optical spectra observed with EOS at HST. 
The high quality of the dataset allowed these authors to perform the 
first physical and chemical analysis of a proplyd, simultaneously 
determining electron densities and temperatures from several indica- 
tors as well as ionic and total abundances of several elements in the 
cusp, tail and jet of LV2. The authors found a high density contrast 
between the cusp (~ 7.9 x 10^ cm"^), the tail (~ 38,000 cm"^) and 
the adjacent nebular background (~ 4700 cm~^), while the tempera- 
tures remain similar for the diff'erent structures. They also found that 
the gaseous heavy-element abundances are generally higher in the 
proplyd compared to the nebular background of the Orion Nebula. 
The positive correlation found between the host star metallicity and 
the presence of giant planetary companions (e.g. [Neves et al.|20"09l 
and references therein) has aroused a great interest in the question 
of the chemical composition of planet-formation circumstellar en- 
velopes. It is therefore important to check if the observed abundance 
difl'erence between LV2 and the Orion Nebula is a general property 
of the proplyds. 

The chemical analysis carried out by [Tsamis "et al. ( 201 1 ) shows 
that a detailed study of proplyd spectra is a way to investigate the 
so-called abundance discrepancy (AD) problem. This is currently 
the most controversial problem in nebular astrophysics (see Liuj 
2002 , Esteban 2002 , Tsamis et al.|2003 , Peimbert & Peimbert 2006 ; 
Liu 2006, Stasinska et al. 2007; Garcia-Rojas & Esteban 2007, 
Rodriguez & Garcia-Rojas^20 10^, Simon-Diaz & Stasihska^2011j ). 



The AD is the disagreement between the abundances of the same ion 
derived from collisionally excited lines (CELs) and recombination 
lines (RLs) of heavy elements. For instance, Garcia-Rojas & Es teban[ 
( [2007| found that, for a large sample of Galactic and extragalactic 
H II regions, the 0^^/H^ ratio calculated from RLs is always between 
0.1 and 0.3 dex higher than the value obtained from CELs for the 
same region. This diff'erence is known as the abundance discrepancy 
factor: ADF(0^^). For the Orion Nebula, numerous studies have 
shown that ADF(0^^) is ~ 0.15 dex for the background nebular gas 
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Figure 1. Section of the central part of the Orion Nebula, a combination of the HST/WFVC2 images taken with different filters ( [O'Dell & Wong 1 1996) . The 
PMAS field of view of 16" x 16" is shown by the white square. The inset zooms in the PMAS image in Ha were rebinned to 160 x 160 pixels and smoothed 
using a Gaussian filter. The red boxes represent the extracted spaxels to analyse the emission from the proplyds 177-341, 170-337 and 170-334, while the yellow 
boxes correspond to their adopted nebular backgrounds. 



(e.g.|Peimbert et al.|1993|[Esteban et al. [19981 [20041 rMesa-Delgadol 
|et al.|2008| l. An initial attempt to investigate the possible role of the 
proplyds in the AD problem fMesa -Delgado et al.|2 00F) did not give 
a clear answer due to the presence of absorption features around 
the O n RLs in the proplyd spectra, the absence of proper high- 
density indicators, and the impossibility of isolating the proplyd 
emission using long-slit spectroscopy (s ee|Henney & O 'Dell 1999}. 
A striking result found by |Tsamis et aL| ( |201 l| t is that the ADF(02+) 



is virtually zero in the cusp due to the high density of the proplyd, 
which produces a clear enhancement of the O^^ abundances from 
CELs with respect to the nebular background, while those from RLs 
are basically similar in both cases. 



The main aim of the present paper is to study several proplyds 
in the Orion Nebula using integral field spectroscopy in the optical 
range at spatial scales of about I" and explore their impact on the 
AD problem. At the same time, we derive the physical conditions 
in the proplyds from empirical diagnostics of temperature, density, 
and ionization. In ^we describe the observations and the data 
reduction procedure. In ^the emission line measurements and the 
reddening correction as well as some representative maps of those 
quantities are presented. In ^we analyse the spatial distributions of 
the excitation ratios, physical conditions and chemical abundances. 
We also attempt to separate the intrinsic emission of the proplyds 
from that of the background nebula. The effect of uncertainties in the 
dust opacity of the proplyd on the background subtraction procedure 
is discussed in ^ In ^we present preliminary results of a physical 
model for the proplyd 177-341 and discuss how these complement 
the more empirical analysis. Finally, in ^our main conclusions are 
summarized. 



2 OBSERVATIONS, DATA REDUCTION AND 
EXTRACTION OF PROPLYD SPECTRA 

2.1 Observations and data reduction 

Integral field spectroscopy was carried out on 2008 December 19 
with the Potsdam Multi- Aperture Spectrometer (PMAS, |Roth et al.| 
|2005| t on the 3.5m telescope at the Calar Alto Observatory (Almeria, 
Spain). The standard lens array integral field unit (IFU) of 16" x 16" 
field of view (FoV) with a spaxel size of x was used. The 
V600 grating was used at two rotator angles to cover the blue part 
of the spectrum from 3500 to 5100 A (angle of -72°) and the red 
part of the spectrum from 5700 to 7200 A (angle of -68°). 

We achieved an eff'ective spectral resolution FWHM ~ 3.6 A in 
both spectral ranges. The total integration times for the blue and red 
spectra were 3 x 400 s and 300 s, respectively. Additional short 10 
second exposures were also taken for each range, avoiding the satu- 
ration of the brightest emission lines. Difi'erent calibration images 
were also taken: HgNe arc lamps for the wavelength calibration and 
continuum lamps needed to extract the 256 individual spectra on the 
CCD. The spectrophotometric standard stars Feige 100, Feige 34, 
Feige 110 and G 191-B2B ( |Oke|1990| l were observed to perform 
the flux calibration. The error of the absolute flux calibration was 
about 5%. The observing night had a few thin cirrus and an average 
seeing of about 1"5. The reduction package specked of iRAi[^was 
used to perform the usual reduction process (bias subtraction, spec- 



^ iRAF is distributed by the National Optical Astronomy Observatory 
(NOAO), which is operated by the Association of Universities for Research in 
Astronomy (AURA) under cooperative agreement with the National Science 
Foundation (NSF). 
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tra extractions and flat-field correction) as well as both wavelength 
and flux calibrations. 

The integral field data are aflfected by the diflferential atmo- 
spheric refraction (DAR), which produces wavelength-dependent 
offsets on the plane of the sky along the local parallactic angle, with 
a magnitude that depends on the air mass of the observation and the 
local atmospheric conditions. We correct for this eff'ect by imple- 
menting the method outlined by Hohenkerk & Sinclair] (19 85 p]in 
our iDL routines, which work in a similar way to those written by 
P.T. Wallace for starlink and jWalsh & Roy| ( [l990| l. The procedure 
calculates fractional spaxel shifts for each wavelength with respect 
to a given reference. We select H/? and Hof as reference wavelengths 
for the blue and red spectral ranges, respectively. The maximum 
DAR shifts for our observations were 0V59 in right ascension and 
V!\5 in declination between [O ii] /13728 and Yip. The emission line 
maps of }Ap and Hof were also used to spatially align the blue and 
red ranges. We found a spaxel oflfset of 0736 in right ascension 
by looking for the maximum correlation between both HyS and Hof 
maps. Taking into account the calculated DAR fractional shifts and 
the alignment oflfset, the PMAS data were finally corrected using 
the BILINEAR function of idl, which uses a bilinear interpolation algo- 
rithm to compute the DAR-corrected data. The final result of this 
process reduces the original FoV to \A" x \A" . The DAR correction 
was tested by checking the alignment of the proplyd peak positions 
at diff'erent wavelengths (e.g., Hi lines), finding a good agreement. 

The observed field is presented in Fig.[T] It is located roughly 
25" to the south east of 6^ Ori C at coordinates: RA = 5^ 35"^ ITA 
and Dec = -5° 23' Al" . The field covers interesting morphological 
structures, in particular several proplyds: 170-334, 170-337, 171- 
340, 173-341, and 177-341, of which 177-341 is clearly the most 
prominent. This field also includes a faint arc of emission in front of 
177-341, in the radial direction to 6^ Ori C. This arc is a bowshock 
produced by the interaction of the photoevaporated flow from the 
proplyd ionization front and the fast wind from the main ionizing 
star (see Henney 2000 and references therein). Another interesting 
feature is the Herbig-Haro object HH 514, which begins as a red- 
shifted microjet arising from the proplyd 170-337. This microjet has 
a counterjet in the opposite direction, but it is too faint to be detected 
in even the HST images (see Bally et al. 2000, O'Dell & Henney] 
^8 ) . In Fig. [T] we have also included the Hof image obtained with 
PMAS and rebinned to 160 x 160 pixels using the rebin function 
of IDL, which performs an expansion of the original map through a 
linear interpolation. We further applied a Gaussian filter to smooth 
the original Hor image. With this method, the spatial distribution of 
certain features observed with the HST in our PMAS maps has been 
nicely reproduced: three out of five proplyds (173-341 and 171-340 
are too faint) and the faint bowshock in front of 177-341. In the 
Hof map a flux enhancement can be also located in the northwest 
corner of the field, which is probably related to the HH 514 object. 
However, part of this information is lost following the DAR cor- 
rection when the FoV is reduced to 14'' x 14". This rebinned and 
smoothed Hor map is used to construct the Hor contours presented 
in all maps throughout this paper. It should be noted that use of 
these contours on the sample maps produces a slight spatial shift 
between the maximum contour value and the maximum intensity 
on the maps. The Hof contours are mainly useful in illustrating the 
localization of the various structures. 



2.2 Extraction of proplyd spectra 

We have extracted several representative spaxels with the aim to 
analyse in detail the nebular properties of the proplyds detected in 
our Hof map. In Fig.[T]the red boxes represent the selected spaxels, 
which coincide with local maxima of Hof emission corresponding 
to the proplyds 177-341, 170-337 and 170-334, located at the co- 
ordinate pos itions (5:4, -2: -i f] (-4 : -5, 5 : 6) and (-5 : -6, 2 : 3), 
respectively. |B ally et al.|(|1998! 2000 ) obtained a wonderful set of 
narrowband images at diff'erent wavelengths ([O iii], [O i], Ha, [Nii] 
and [Sii]) for several proplyds, including 170-337 and 177-341. 
These images show that the maximum emission line intensity comes 
from their bright cusps, while the tails are much fainter. Given the 
very small angular sizes of the proplyds, it is clear that the observed 
spectra of the areas of the field we associate with proplyds are only 
a composite of both proplyd and background emissions. The largest 
proplyd is 177-341, with a diameter of about 0'.7, but even this is 
smaller than our spatial resolution and the average seeing during 
the observations. Therefore, it is vital to disentangle the intrinsic 
emission of these objects from that of the ambient ionized gas and 
other possi ble contributions ([Henney & 0'Dell|1999[|Vasconcelos| 
|et al. 12005) [Tsamis et al.|201 l |l. To do that, we have selected a single 
spaxel as representative for each proplyd, ensuring that it contains 
the entire or most of the emission from their cusps. In principle, the 
main source of contamination is the emission from the background 
of the nebula, which has to be properly subtracted. However, it is 
diflficult to define an accurate nebular background since the neb- 
ula shows significant brightness variations at small spatial scales. 
Furthermore, the opacity due to the dust inside the proplyds is an- 
other unknown factor that may introduce systematic errors in the 
subtraction process of the nebular emission (see This issue is 
extensively discussed in [Henney & 0'DeITHr999 ), who found an 
eff'ective extinction of Ay ^ 0.1-0.2 in the proplyd 177-341, which 
implies eff'ective extinction coefficients for the proplyd, c(H5)p, of 
about 0.04-0.08 dex for their observing conditions and extracted 
aperture. 

In order to put limits on the possible effects of proplyd dust on 
our extracted spectra, we consider two extreme possibilities. First, 
the case where the proplyd lies in front of all the nebular emission 
and is totally opaque, c(H/?)p » 1, so that the observed (background 
plus proplyd) emission directly corresponds with the proplyd emis- 
sion itself. Second, the case where the proplyd is totally transparent, 
c(I^)P = 0, or lies behind the nebular emission, so that the proplyd 
spectrum must be obtained by subtracting the full nebular emission 
from the observed spectrum. These two cases are particularly easy 
to deal with since they do not depend on the properties of the dust in- 
side the proplyds. The more complex intermediate case of a partially 
translucent proplyd is deferred until ^ 

Taking into account the above considerations, it seems rea- 
sonable to obtain the intrinsic proplyd spectra by assuming the 
transparent case, that is a direct subtraction of the nebular back- 
ground from the observed spectra (proplyd plus background). This 
is the same approach that was applied in integral field studies of 
the LV2 proplyd ( [Vasconcelos et al.|2005| [Tsamis et al.|2011[ l. As 
previously mentioned, the main diflftculty with this method lies in 
defining the nebular background. As a result, only rough estimations 
can be obtained from observations assuming an homogeneous back- 
ground around the proplyds. In this case, the background emission 



^ Reproduced, with permission, from data supplied by KM Nautical Al- 
manac Office, UKHO © Crown Copyright. 



^ Hereinafter, the localization of an individual spaxel or multiple spaxels 
will be given as {a\ :a2,6\ :52) according to the coordinate system used in 
the maps shown in next sections. 
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Figure 2. The PMAS spectra of the proplyd 177-341 [spaxel position 
(5 : 4, -2 : -1)] from 4600 to 4700 A, showing the On Ml multiplet lines. 
The black (upper) line shows the observed spectrum (proplyd plus back- 
ground emission), while the blue (middle) line is the average spectrum of 
the background emission (see Fig.[^. The red (lower) line is obtained as a 
direct subtraction of the average background emission from the observed 
one (see ^2.2] representing the spectrum of the proplyd 177-341 itself. The 
flux of the proplyd has been multiplied by 4 to facilitate the comparison with 
the other spectra. 

was obtained as an average spectra from several spaxels around each 
proplyd. These spaxels are represented as yellow boxes in Fig. [T] 
The proplyds 170-334 and 170-337 shared the three spaxels to the 
east of 170-337 to construct their associated background spectra. In 
order to check the consistency of our results, additional background 
spectra were obtained from other parts of the field to compare them 
with the backgrounds assigned to the proplyds. We estimate that 
diflferences among the diflferent backgrounds are lower than 10%. 
This variation in the background intensity contributes significant 
uncertainty to the proplyd flux measurements of those lines that are 
either (1) relatively weak in the proplyd compared with the back- 
ground (e.g., [S ii] nebular lines; see ^3.1| l, or (2) relatively faint 
compared with the adjacent continuum emission (e.g, Cii or On 
RLs). 

Finally, we also have to consider other possible sources of 
contamination in the spectra of the selected spaxels, such as the 
presence of high- velocity emission from any jet associated with the 
proplyds. Only for 170-337 has a jet has been reported, HH 514 
( [Bally et al.|1998l[2000l|OT)ell & Henn ey 2008 ), but it is rather faint 
and its contribution to the total emission should not be important. 



3 LINE MEASUREMENTS AND REDDENING 
CORRECTION 

3.1 Line measurements 

The flux of each emission line for each individual spaxel in the 
PMAS field was measured by fitting a Gaussian profile plus a linear 
continuum using the splot routine of iraf. Where necessary, such as 
for blended lines, multiple Gaussian components were fit simultane- 
ously. The brighter lines were measured by an automated process 
using our own iraf scripts, while the fainter ones (e.g., C ii A4261 
or O II /14650) were measured by hand in all spaxels. The emission 
line fluxes of the individual proplyd spectra (background- subtracted 
and observed) and of their associated backgrounds were also mea- 
sured by hand using the same methodology. The errors in the flux 



measurements were determined according to the criteria defined 
by Mesa-Delgado et al. (2008 ). The total error for each line was 
estimated as the quadratic sum of the flux measurement error and 
the flux calibration error. 

In the red spectral range we found a disagreement between 
the emission fluxes measured in the short and long exposures. It 
is important to correct for this eflfect because the brightest lines 
of that range, such as Ha, are saturated in the long exposures and 
must be are measured in the short exposures. The disagreement is 
insignificant (lower than the calibration error) in most spaxels cov- 
ering the nebular background, but is larger at the proplyd positions. 
The measured short-exposure fluxes are on average 11%, 8% and 
3% higher than the long-exposure fluxes for 177-341, 170-337 and 
170-334, respectively, which leads to an even larger discrepancy in 
the background- subtracted values of 35%, 29% and 18%, respec- 
tively. Since this only afl'ects the point sources within the FoV, it is 
likely due to seeing variations during the observing night. Indeed, 
the eff'ect largely disappears (becoming smaller than the calibration 
error) if a larger sample region is used to obtain the intrinsic proplyd 
emission of 177-341. We therefore correct for the eff'ect by using the 
above factors, which are rather constant throughout the red spectral 
range, to rescale the short-exposure fluxes. For the blue spectral 
range the diflferences between the short and long exposures were 
always lower than the calibration error, and so no correction was 
necessary. 

The emission lines used in our analysis were: Hi Balmer lines, 
from Ha to HI 2, to compute the reddening correction and test the 
DAR correction; CELs of difl'erent ionic species to derive physical 
conditions and chemical abundances; and the faint C ii and O ii lines 
to determine the ionic abundances from RLs. These lines were easily 
detected in the whole field, excluding the extracted spaxels of the 
proplyds 170-337 and 170-334, and a few spaxels around them. The 
spectrum of the proplyd 177-341 was the only one with high enough 
signal-to-noise ratio to detect these faint lines in both observed and 
background-subtracted spectra. A section of the extracted spectra of 
this proplyd is presented in Fig. [2] where we can see the blend of 
O II /14649, /14651 lines before and after the background subtraction. 

In Fig.[3]we present several representative line flux maps. The 
diflferent rows show the spatial distributions of auroral (left-hand 
column) and nebular (right-hand column) emission lines for the 
same ion: S+ in Figs[3];a)-[3];b), in Figs [3];c)-[3];d) and 0^+ in 
Figs[3|e)-[3jf). For all three proplyds, there are significant diflfer- 
ences between the spatial distributions of the auroral and nebular 
lines. The [S ii] flux distribution shows the most extreme behaviour, 
whereby the proplyds are barely detectable in the [S ii] A61\l nebu- 
lar line map, except perhaps for 170-337 (Fig.jSj)). In contrast, they 
show clear local maxima in the distribution of [S ii] /14068 auroral 
line (Fig. [3^). In fact, for 177-341, we do not detect any trace of 
[Sii] A6711 emission in the background- subtracted spectrum. On 
the other hand, the spatial distributions of [N ii] and [O iii] lines 
show a somewhat diff'erent behaviour for the diff'erent proplyds. The 
[N ii] A5155 map has prominent peaks at the positions of the three 
proplyds (Fig.|3j:), whereas the [Nii] nebular emission (Fig.[3ji) is 
more extended at the positions of the proplyds 170-334 and 170-337, 
and has a secondary maximum in 177-341. For [Oiii], the auroral 
line map (Fig. [3^) has its maximum at the position of 170-334, a 
secondary peak at 177-341 and an extended distribution at 170-337. 
The spatial distribution of the [O iii] nebular lines (Fig.|3}i) seems to 
be somewhat more extended at the positions of the three proplyds. 
We have found similar behaviour for other emission lines. For in- 
stance, [O i] /16300 has a spatial distribution very similar to [N ii] 
/15755 and [S ii] /14068 with prominent peaks at the positions of 
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Figure 3. Spatial distribution maps of auroral (left-hand column) and nebular (right-hand column) emission lines of different ions: a, b) [S ii] /14068 and A61\l\ 
c, d) [Nil] A5155 and /16583; and e, f) [Oiii] /14363 and AA959. The fluxes are in units of erg cm~^ s~^ and Hof contours are overplotted in all maps. 



the proplyds. Another interesting case is the [O ii] /13728 emission 
line, whose spatial distribution is very similar to that of the [S ii] 
nebular line, showing the same absence of peaks at the proplyd 
locations. Indeed, we did not detect any emission of [On] /13728 
line at 177-341 after the background subtraction. Other emission 
lines such as [Neiii], [Siii] and [Arm] present spatial distributions 
similar to those of the [Oiii] and [Nii] nebular lines. 

In Table [T] for each of the main CELs we show the ratio of the 
background- subtracted proplyd flux to the total observed flux in the 
extracted spaxels (see |2.2| >, with the aim of identifying from which 
component nebular and auroral lines are emitted. Regarding the 
emission lines presented in Fig.|3] we see that the [S ii] auroral lines 
in the background- subtracted spectra of the proplyd 177-341 repre- 
sent about 75% of the observed flux, as opposed to only 5% for the 
nebular lines. For [Nii], the nebular lines are dominated by the back- 
ground emission (76-88% of observed flux), whereas the proplyd 
emission dominates (56-77%) in the auroral line. Combining the 
information of Fig.[3]and Table[T] it is found that, compared to the 
Orion background, 177-341 has a much higher surface brightness in 
the [S ii] /14068 auroral line than in the [S ii] X61\l nebular line (fac- 
tor of 15), followed by 170-334 (factor of 9.8) and 170-337 (factor 
of 3.2). A similar pattern is seen for both [Nii] auroral and nebular 
lines (with more similar factors of 3.2, 2.5 and 4.7, respectively, for 
the three proplyds). In contrast to the [S ii] and [Nii] lines, where 
the emission of nebular (auroral) lines from the proplyds is a small 
(large) fraction of the total observed (proplyd plus background) 




6 4 2 -2-4-6 
Aa (arcseconds) 

Figure 4. Spatial distribution map of the extinction coefficient, c(l^), along 
the Une of sight with Ha contours overplotted. 



emission along the line of sight, the [O iii] auroral and nebular line 
fraction about 40-50% and 25% respectively. As we will discuss in 
g4.2.1| these patterns can be attributed to collisional de-excitation at 
the high electron densities expected in the proplyds (up to 10^ cm~^; 
[Henney et al.|2002|[Tsamis et al.|20lT] ). 
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Table 1. Critical densities^, yf^^, of the listed emission lines and the fraction (in %) 
of the line fluxes that come from the proplyd itself (background- subtracted spectrum) 
with respect to the observed spectrum (proplyd plus background emission). 



F(proplyd)/F(observed) (%) 



Ion 


X{k) 


nf' (cm-3) 


177-341 


170-337 


170-334 


[Nn] 


6548, 6584 


9x10"^ 


24 


24 


12 




5755 


10^ 


77 


61 


56 


[On] 


3726, 3729 


10^,4x103 





17 


9 




7318, 7329 


4x10^ 








[Om] 


4959, 5007 


7x10^ 


24 


22 


27 




4363 


2x10^ 


37 


38 


51 


[Neni] 


3869 


10^ 


17 


35 


27 


[Sn] 


6717, 6731 


2x103,4x10^ 


5 


17 


4 




4068, 4076 


8x10^2x10^ 


75 


54 


39 


[Sm] 


6312 


10^ 


41 


32 


32 


[Arm] 


7135 


5x10^ 


28 


30 


26 



For the atomic data of |Garcia-Rojas et al.H2009) . 



3.2 Reddening correction 

The extinction coefficient, c(HyS), along the line of sight was derived 
for each spectrum by comparing the observed YiyfiAp and Yid/Yip 
ratios with the theoretical values predicted by [Storey & Hummer] 
( |1995| | for = 1000 cm-^ and = 10,000 K. Additional H i line 
ratios were also measured but they were discarded in deriving the 
final c(Rp) value because they were rather faint (Hll and HI 2) 
or were blended with other brighter lines (H7, H8 and H9). For 
the background- subtracted proplyd spectra we used the same ratios 
despite the very high densities expected for the proplyds. We do 
not expect significant changes in the results because of the weak 
density-dependence of the Balmer line ratios. In fact, the selected 
ratios show a variation less than 2% when electron density and tem- 
perature change from 1000 to 10^ cm-^ and from 7500 to 12,500 K 
( [Storey & Hu mmer 1995 |. The extin ction law normalized to H/?, 
/(/I), determined by Blagrave e t al.| ( [2QQ7| l for the Orion Nebula 
{Ry - 5.5) was adopted. This extinction function produces c(H/?) 
coefficients about 0.1 dex higher than the classical one derived by 
( [Costero & Peimbert|1970j l. For each spectrum the final c(H/?) value 
was derived as an average of the individual c(H/?) determinations 
from each Balmer line ratio, weighted by their uncertainties, giving 
results with a typical uncertainty in c(H/5) of about 0.20 dex. All flux 
measurements were referenced to H/? and Hof for the blue and red 
range, respectively. After the reddening correction, we produce a fi- 
nal homogeneous set of flux ratios by re-scaling the flux ratios of the 
red range to H/? using the theoretical Hof/H/? ratio for - 10,000 K 
and Hq - 1000 cm"^. In Table|2]we present the dereddened flux ra- 
tios of the main emission lines used in our analysis for the extracted 
proplyd spectra together with their associated backgrounds. 

The spatial distribution map of c(HyS) along the line of sight is 
shown in Fig.|4] The c(H/?) coefficient ranges from 0.5 to 1.0 dex, 
with an average value of 0.74 dex and a standard deviation of 0.14 
dex. This mean value is consistent with the measurements of O'Dell 
[& Yusef-Zadehl ( [2000l t, who derived the c(H/?) map from the Hof/H/? 
line ratio, finding values of 0.6 to 0.8 at the position of our PMAS 
field, even though a diflferent extinction law was used by these au- 
thors. In Tables[3]and|4]we present the extinction coefficients along 
the line of sight, labeled as c(H/?)^°% obtained for the background- 
subtracted spectra (c(I^)P = 0), the observed spectra (c(H/?)p » 1) 
and the adjacent background for each extracted proplyd. From the 
extinction map we have estimated mean c(H/?) coefficients at the 
proplyd positions, which amount to 0.80 dex for 177-341, 0.83 dex 
for 170-337 and 0.80 dex for 170-334. We find a good agreement 



for each proplyd between these values and those of Tables |3] and 
|4| considering the uncertainties. The background values also agree 
with the extinction map of Fig.|4| Furthermore, our extinction values 
along the line of sight of the proplyds 177-341 and 170-337 are also 
consistent with those presented by |0'Dell[fl998| t, 0.81 dex and 0.84 
dex, respectively. Finally, we note that the reddening determinations 
for the background-subtracted and observed spectra of each proplyd 
are in agreement considering the errors. This is expected if both 
emissions arise from the same spatial positions along the line of 
sight. However, the small diflferences between those pairs of values 
could be the result of internal extinction produced by the presence 
of dust inside the nebula or the contribution of scattered light. 



4 RESULTS AND DISCUSSION 

4.1 Spatial distributions of line flux ratios 

In this section we explore the behaviour of the [NiiJ/Hof and 
[O ni]/H/? line ratios to study the spatial distribution of the gas exci- 
tation. Assuming that the pattern observed in the spatial distributions 
of Fig.[3]is produced by the effect of collisional de-excitations (see 
^4.2. lj >, we have plotted in Fig. [5] the excitation ratios obtained from 
nebular and auroral lines. On the one hand, we show the typical 
ratios using the [N ii] i6548, i6583 and [O m] ^4959, .15007 nebular 
lines (Figs [5^ and[5]:), which allow us to analyse the excitation of 
the background emission. On the other hand, we have also used 
the [Nil] .15755 and [Om] .14363 auroral lines (Figs|5]D and[5}i) to 
explore the excitation ratios at the proplyd positions, where the neb- 
ular lines would be aff'ected by collisional de-excitation efi'ects. The 
ionization structure of the background emission seems to be rather 
similar using both nebular and auroral line ratios. The low excitation 
gas is basically surrounding the high excitation areas as we see in 
Figs[5|a) and[5jc). The maximum emission of the low ionization 
gas comes from the northwest corner, while the high ionization gas 
is located at the south side of the field. This stratification seems to 
be related to the edge of the Red Bay and some of the Dark Lanes, 
which are kinematic features identified bylGarcia-Diaz & H enneyj 
( [2007| in this area of the Orion Nebula. In the case of high ionization 
gas, some of the higher values seem to be located at the position 
of the bowshock, spaxels (2 : 0, -3 : 0). At the proplyd positions, 
the excitation ratios using the nebular lines reveal that the emission 
from the proplyds 177-341 and 170-337 comes from low ionization 
gas. This is clearly confirmed when we use the excitation ratios 
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based on auroral lines, Figs[5|b) and|5|d), where we find prominent 
spikes at these positions. It is an expected behaviour because we are 
observing the radiation emitted from the low ionization gas in the 
ionization front of the cocoon structure assumed for the proplyds 
(see figure 5 of |Henney & 0'Dell|1999] . The proplyd 170-334 also 
has low ionization emission, but the emission of auroral lines of high 
ionization gas is dominant for this proplyd at spaxels (-4 : -3, 5 : 6). 
This behaviour may imply a extremely high density for the pro- 



plyd 170-334, which is quenching the emission of the [Om] /14959, 
A5001 nebular lines (see 



4.2 Physical conditions 

In this section we present electron densities, n^, and temperatures, 
Te, derived from the analysis of each spaxel as well as the extracted 
spectra of the proplyds. The spatial distribution maps of and 
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Uq are shown in Fig. [6| and the physical conditions derived for the 
background- subtracted, observed and adjacent background spectra 
of each proplyd are available in Tables [3] and |4| In Table [3] we 
also present the different results of the proplyd 177-341 by taking 
into account the effect of opacity in the proplyd on the background 
subtraction procedure (see Here the opacity is parametrized by 

The calculations for each spaxel were performed using the 
TEMDEN task of the nebular package ( |Shaw & Dufour||1995|l of 
iRAF with updated atomic data (see table 6 of Garcia-Rojas et al.] 
|2009| l. An iterative procedure was used to determine and 
from the usual CEL diagnostic ratios: [S ii] /16717//16731 for n^, and 
[Oiii] ^^4959,5007M4363 and [Nii] /1/16548, 6583M5755 for T^. 
We assume an initial - 10,000 K to derive a first estimate foe 
fZe([S ii]), then re([0 iii]) and re([Nii]) are calculated. The process is 
repeated until convergence, which is reached after 3 or 5 iterations, 
taking re([Nii]) as representative for ion assuming a two-zone 
scheme (see |4.3| l. 

The physical conditions given in Tables |3] and |4] for the ex- 
tracted spectra of the proplyds (background- subtracted and ob- 
served) as well as their adjacent backgrounds were computed fol- 
lowing a similar iterative procedure. Additionally, we use alternative 
density diagnostics that are more suitable for deriving the expected 
high electron densiti es for the proplyds, which can reach values 
exceeding 10^ cm"^ ( [Henney et al.|2002||Tsamis et al.|201l| . At 
those densities, diagnostics based on typical optical line ratios (e.g., 
[Sii] ^6717M6731 or [Cliii] ^5517M5537) would not return reU- 
able values due to the low critical densities of the upper levels of 
the transitions producing those lines. The first alternative density 
indicator is based on the analysis of [Feiii] emission lines, which 
are useful as robust density diagnostics with much wider ranges of 
validity than other line ratios. As shown in Table[2] we have detected 
several [Fe iii] lines (a total number of 5 lines of the 3F multiplet 
and one of the 2F multiplet) in the background- subtracted spectra 
of 177-341, in the observed spectra of 177-341 and 170-337, and in 



the three adjacent background spectra of each proplyd. The density 
adopted from the analysis of [Feiii] lines, ;ie([Feiii]), is the value 
that provides the minimum dispersion between the observed and 
theoretical line ratios with respect to [Feiii] /14658 for the whole set 
of lines. The theoretical [Fe iii] emissivities were computed by solv- 
ing a 34-level atomic model using the atomic data of|Zhang|( 1996| 
for collisional strengths, and ' Quinet| (|1996^ and ' Johansson et al.| 
(2000 ) for transition probabilities. The second alternative density 
indicator consists in obtaining diagnostic curves in the plane Uq-Tq 
from the line ratios [Nii] ^6583/^5755, [Oiii] ^5007/^4363, [Sii] 
^^6717,6731/^^4068,4076 and [Sii] ^6717/^6731. The se diagnos- 
tic curves (Fig.[7| were computed making use of pyneb ( Luri diana| 
letaLpOTTT), which is an updated version of the nebular package 
written in python programming language, and the set of atomic data 
presented in ,Garcia-Rojas et al.,^2009j . In general, the intersection 
of those curves yields a clear solution for Uq-Tq, whose values are 
labeled as diag. Uq-T^ in Tables |3] and |4] In particular, we have 
chosen the values provided by the [Nii] and [O iii] curves. The [S ii] 
curves obtained for the background- subtracted and observed spectra 
of the proplyds require a more detailed discussion because apart 
from collisional de-excitation of the [S ii] nebular emission lines, 
uncertainties in the background subtraction could be affecting the 
results (see next section). Finally, in the cases where both previous 
density indicators were available for the background-subtracted and 
observed proplyd spectra, the adopted values were calculated as 
weighted averages, otherwise the result from the diagnostic curves 
was adopted. 

For the adjacent background spectra it was not necessary to use 
the diagnostic diagrams and the densities derived from the [Fe iii] 
and [S ii] lines were weighted to obtain the average final values. 
Then, re([Oiii]) and re([Nii]) were calculated using these adopted 
densities. We did not correct for the recombination contribution 
to the flux of the [Nii] /15755 auroral line in any of the above 
calculations. For the background emission and the observed spectra 
this is expected to be rather low because of the ionization degree of 
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Figure 6. Spatial distribution maps of physical conditions with Hof contours overplotted: a) /ie([S n]), b) re([Nii]), c) re([Oin]) and d) re([Nii])/re([Oiii]) ratio. 



the Orion Nebula (e.gJ Esteban et al.|2004| l and even less important 
in the case of the proplyds. 



4.2.1 Electron densities 

In the spatial distribution map of ^e([S n]) shown in Fig.[6ja), we do 
not find sharp density contrast between the proplyds and background. 
The typical values range from 3000 to 8000 cm~^ with uncertainties 
between 1000 and 2000 cm~^. On average the nebular background 
has a density of about 3700 cm~^ ranging from 2800 to 5000 cm~^. 
The maximum value of We([S ii]) (- 8000 cm~^) is reached at the 
position of the proplyd 170-334 and secondary peaks of about 5800 
and 5200 cm'^ are found at the proplyds 170-337 and 177-341, 
respectively. As is expected, these values are in agreement within 
the uncertainties with «e([S ii]) values shown in Tables |3]and|4]for 
the observed extracted spectra of the proplyds: 3700 cm~^ for 177- 
341, 4240 cm-3 for 170-337 and 8540 cm-^ for 170-334. However, 
the density determinations obtained from [S ii] line ratios clearly 
disagree with those derived from the [Fe m] lines and the diagnostic 
curves also included in Tables [3] and |4] The two last indicators point 
out the existence of high-density ionized gas at the proplyd positions, 
where the values derived for the observed spectra (proplyd plus 
background) now amount to 93,100 cm~^ for 177-341, 57,900 cm~^ 
for 170-337 and 24,000 cm-^ for 170-334. 

We interpret this disagreement as due to the fact that the [S n] 
/16717, /16731 emission detected in the observed spectra of the pro- 
plyds comes mainly from the nebular background. We saw in Table[T] 
that the proplyd emission itself only represents 4-17% of the total 
observed emission of the [S ii] nebular lines. If the high densities 
derived from the [Fe iii] lines and the diagnostic curves correctly 
reflect the proplyd physical conditions, then the [S n] line ratio will 
not provide reliable densities because collisional de-excitation will 
be severely suppressing the emission of [S ii] nebular lines emitted 
from these objects due to the low critical densities of the upper lev- 
els (see Table [T}. Therefore, the ^e([S n]) values determined for the 
observed spectra are more similar to the background ones and this 



can be seen in Tables [3]and[4| This also explains the behaviour of the 
[S ii] /16717//16731 curves of the observed spectra of the proplyds 
(Figs[7]3,|7^ and Th), which are similar to those of the background 
emission (Figs|7 and|7}). Another explanation for this would be 
that the [S ii] nebular lines could partly originate from a lower den- 
sity envelope surrounding the proplyds, but our observations do not 
have enough spatial resolution to confirm this; however, the generic 
Orion emission is most probably the dominant contributor to these 
lines (see Table [TJ. Then, the discussion above clearly indicates that 
the [S ii] AA6111, 6731//1/14068, 4076 diagnostic curves are unsuit- 
able for deriving densities from the observed spectra of the proplyds 
because the constituent lines of the aforementioned [S ii] nebular 
to auroral diagnostic ratio arise in two very diflPerent gas volumes 
(generic nebula vs. proplyd). Thus in this case the ratio has a very 
limited diagnostic value. In the same way as for the [S ii] nebular 
lines, collisional de-excitation will aflfect very strongly the flux of 
the [O ii] /13728 doublet due to its similarly low critical densities 
(see Table [T}. 

The situation discussed in the previous paragraph hampers the 
correct interpretation of the background- subtracted spectra. The 
[Sii] /16717, /16731 and [On] /13728 emission should partially or 
totally disappear after subtracting the background emission from 
the observed spectra. However, these lines do not completely dis- 
appear in the background-subtracted spectra of 170-337 as we see 
in Table |2] In fact, the relative contribution of the proplyd to the 
[S ii] and [O ii] emission is larger for this object than for the others 
(see Table [TJ. Furthermore, we derive a density of about 2300 cm~^ 
for 170-337 using the [S ii] line ratio, which is actually more con- 
sistent with the background values and very diflferent from the high 
density determined from the diagnostic curves, 1.29 x 10^ cm~^. 
This is also the behaviour that we observe in the diagnostic dia- 
gram obtained for the background- subtracted spectra of this proplyd 
(Fig. |7ji), where the [S ii] curves provide lower densities than the 
intersecion of [Nii] and [Oiii]. These results indicate that there is 
residual background contamination in the background- subtracted 
spectra of this proplyd, probably due to its small size. We expect a 
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Figure 7. Uq-Tq diagnostic curves for the background-subtracted, observed and adjacent background spectra for the extracted proplyds: 177-341 (first row), 
170-337 (second row) and 170-334 (third row). The different curves represent the following diagnostic ratios: [N ii] /16583//15755, [O iii] /15007//14363, [S ii] 
AA6111, 6731//1/14068, 4076 (narrow band) and [S n] /16717//16731 (wide band). The width of each curve represents the diagnostic uncertainty. 



similar situation for the proplyd 170-334, which is even smaller than 
170-337. We do not detect [S n] A61\l in the background-subtracted 
spectrum of this object, but there is some emission from [O ii] /13728 
(see Table [2| and the curves of [S ii], [Nii] and [Om] do not have 
a common intersection (Fig. |7^). Hence, the physical conditions 
of these proplyds should be considered with some caution because 
they could be contaminated by background emission. On the other 
hand, background contamination seems to not be substantial in the 
background- subtracted spectra of 177-341 due to its larger size, 
the optimal signal-to-noise relation of its spectra, and the appro- 
priate isolation of the selected spaxel. The plasma diagnostic of 
this proplyd (Fig. [7^) produces consistent physical conditions for 
the different curves. All these considerations are relevant for the 
analysis of the proplyd properties and, in particular, they have to 
be considered in the calculations of the chemical abundances from 
CELs, as we will see in g4.3| 

It is also interesting to compare the densities of the background 
and proplyds derived from the [Fe ni] analysis and the diagnostic 
curves (Tables [3] and |4]l. We see that proplyd densities are one to 
two orders of magnitude higher than the background ones, a similar 
result to that found for the LV2 proplyd (Tsamis et al. 201 1 ). For 
instance, the proplyd 177-341 has densities of about 3.8 x 10^ cm~^ 
and 9.3 x 10^ cm~^ determined from the background- subtracted 
and observed spectra, respectively, whilst the background density 
amounts to 3800 cm~^. For the other two objects, 170-337 and 170- 
334, we find a similar contrast, although the proplyd densities are 
lower than those of 177-341, perhaps because they are more affected 



by background contamination as we mentioned in the previous 
paragraph. 

In their study of the LV2 proplyd, [Tsamis & Walshl pOTT) 
found a clear correlation between the critical density and the relative 
contribution of the proplyd emission to the total observed one in 
the area occupied by the object. In our case, such correlation can be 
seen in Table [T] although it is not so clear as in LV2. However, the 
densities of our proplyds are different and the background contami- 
nation can also affect the smaller proplyds. The general behaviour 
of the proplyds observed in our PMAS field shows that emission 
lines emitted by ions with high ionization potentials (O^^, Ne^^, 
S^^ and Ar^^) are dominated by the background emission, whereas 
the proplyd contribution to emission lines arising from ions with 
low ionization potentials (N^, and O^) varies depending on the 
transition type: background emission dominates the nebular lines 
and proplyd emission dominates the auroral lines. 

Finally, inspection of Table [T] also gives further insights into 
the effect of collisional de-excitation due to high densities in the 
different spaxels containing proplyd emission. In all cases, we see 
a very strong or even total suppression of the [S n], [O ii] and [Nii] 
nebular lines. This indicates that the true electron densities must 
be above the critical density of the upper levels that produce those 
lines. However, 170-334 presents a somewhat different collisional 
de-excitation pattern because the contribution of the proplyd to the 
flux of [Om] auroral line is higher (see Table [T] and Fig. [3^). In 
addition, the spatial distribution maps of [O ni]/H/? ratio (Fig. |5ji) 
and re([Oin]) (Fig. [6]:) show a clear increase with respect to the 
surroundings at the position of this proplyd. These two facts indicate 
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that the true electron density in this proplyd should be even higher 
than the critical density of the nebular [O ni] lines, i.e., higher than 
7 X 10^ cm~^, much larger than the value quoted for the proplyd in 



Tableg 



of about 1.6 X 10 cm ^. Indeed, the models of Henney & 



[Arthur (1998]) give a peak density for 170-334 (erroneously named 
171-334 in their Table 2) of about 3 x 10^ cm~^, which is 5 times 
higher than that predicted for 177-341 (6.2 x 10^ cm~^). This is con- 
sistent with the [O ni] nebular lines being collisionally de-excited in 
170-334. Probably, the aforementioned background contamination 
is producing the artificial low density value presented in Table |4] 

4.2.2 Electron temperatures 

If we know the true representative of the ionized gas in the pro- 
ply ds, then we can derive the appropriate because statistical 
equilibrium equations take into account the collisional de-excitation 
eflfects. However, this is not so simple because the physical condi- 
tions in a proplyd could be far from homogeneous and, in particular, 
the density is expected to show significative structure at spatial 
scales much smaller than can be resolved with our observations. 
Evidence of such eff'ects was provided by |Tsamis & Walsh] ( [20 ll| l 
who combined high-spectral resolution VLT observations of optical 
[O ni] nebular and auroral lines with the HST far-UV C ni] obser- 
vations by [Henney^ et al.[ ( |2002) l to simultaneously determine the 
electron temperature and density in the O^^ zone of LV2: the tem- 
perature shows a maximum root-mean- square variation of 10% and 
the density varies by up to 50% across the line profile associated 
with the kinematic core of the proplyd. 

Assuming that the adopted densities from [Feiii] lines and 
the diagnostic curves are the correct, we determine re([Nii]) and 
re([0 m]) for both background-subtracted and observed spectra of 
the proplyds (Tables [3]and[4]) and compare these results with the spa- 
tial distribution maps calculated using the ^e([S n]) map of Fig.[6ja). 
The background shows almost constant temperatures, 9800 K for 
re([Nn]) and 8650 K for re([Oni]), with typical errors of about 
600 and 200 K, respectively. In both temperature distribution maps 
(Figs[6j) and[6]:) the proplyds show very prominent spikes, which 
are more evident in the re([Nn]) distribution in the cases of 177-341 
and 170-337 and in re([Oni]) for 170-334. It is clear that these 
peaks do not correspond to the true temperatures of the proplyds 
because the density used to derive the temperature maps, /ie([S n]), 
is not reliable for the proplyds as we previously commented ( |4.2.1| |. 
In Tables [3] and |4| where we have used the higher densities given by 
the other indicators, we can see that the resultant temperatures do 
not show such high contrast. 

For 177-341, re([Oni]) and re([Nii]) are lower than the back- 
ground temperatures by only a few hundred K and about 1000 K, 
respectively. Interpretation of the temperatures in the proplyds 170- 
337 and 170-334 is more diflficult because of the background sub- 
traction problems as we previously mentioned (see ^4.2.1[ ). High 
temperatures derived from the background- subtracted spectra of 
these proplyds are probably due to errors in correcting for colli- 
sional de-excitation eff"ects occasioned by the adopted density de- 
rived from the [Fe iii] lines and diagnostic diagrams being lower 
than the true value. This eff'ect is clear for the proplyd 170-334, 
where the temperatures are even higher than those of the nebular 
background. 

It is likely that much of the apparent structure seen in tempera- 
ture indicator maps, especially re([Nii]) (Fig.[6|3), is in reality due 
to high-density ionized gas. In the spatial distribution of re([Nii]) 
we find two secondary maxima of about 1 1,500 K at the spaxel 
positions (0 : -2, -3 : -2) and (-4 : -5, -2 : 0). These positions do not 



correspond to any noticeable morphological structures in our emis- 
sion line maps (see Fig. [3]), but they do coincide with the positions 
of the faint proplyds 173-341 and 171-340 (see the HST image of 
Fig.[T]). Careful examination of the spatial distribution of the [Nii] 
/15755 line (Fig. [5}:) reveals a slight increase in the line flux at these 
spaxels of about 13% and 20%, respectively, with respect to the 
adjacent background. On the other hand, no similar enhancements 
were detected in the spatial distribution of the [N ii] /16548 nebular 
line (Fig.[3ji). Therefore, re([Nii]) is again overestimated in our 
detected proplyds because the efi'ect of collisional de-excitation is 
not well corrected for due to the adoption of the density derived 
from the [S ii] line ratio. 



4.3 Chemical abundances from CELs 



Making use of the ionic task of the nebular package ( Shaw & D ufour 
|1995|l and the set of atomic data presented in Table 6 of Garcia-Rojas 



[et al.[ (|2009 ), we have derived ionic abundances for several ions from 
CELs: N+, 0+, 0^+, S+, 8^+, Ne^^ and Ar^^. We have assumed a 
two-zone scheme with the absence of temperature fluctuations, t^= 
0, adopting re([Nii]) for ions with low ionization potentials (N^, 
and S^) and re([Oiii]) for ions with high ionization potentials (O^^, 
S^^, Ne^^ and Ar^^). We have also determined the Fe^^ abundances 
for the extracted spaxels using re([Nn]) and the atomic model 
described in ^4.2[ The abundance uncertainties are the quadratic 
sum of independent contributions from n^, Tq, and line flux errors. 

In Fig. [8] we present the spatial distribution maps of the abun- 
dances of two relevant ions, and O^^, as well as of the total O 
abundance calculated as the sum of and O^^ abundances. The 
typical uncertainties for each spaxel are about 0.15 dex for 0^/H+, 
0.04 dex for O^+ZH^ and 0.10 dex for 0/H. In all maps, the mini- 
mum values are clearly located at the proplyd positions and their 
surroundings. Considering the results discussed in ^4.2| this must 
be related to the overestimation of temperatures due to the use of 
an unreliable density for the proplyds. Excluding the proplyd emis- 
sion areas, the mean value for the total O abundance amounts to 
8.45 + 0.05 dex, which is in agreement with the previous determina- 
tions for the Orion Nebula available in the literature (8.50 dex; e.g. 
[Esteban et al.|2004| l. In the spatial distribution of the 0+/H+ ratio 
(Fig. |8^), the values of about 7.6 dex slightly to the south of the 
center of the field correspond to the areas of higher re([Nii]) asso- 
ciated with the proplyd 173-341 at spaxels (0: -2, -3 : -2) and the 
proplyd 171-340 at (-4: -5, -2:0). In Fig.[^c), the total 0/H ratio 
also shows secondary minima at those positions where we found the 
re([Nii]) and re([0 iii]) peaks for the proplyds. The particular struc- 
ture of the high-excitation bowshock in front of 177-341, spaxels 
(2:0,-3:0), does not show any noticeable behaviour in the oxygen 
abundance maps compared with its surroundings. We measure a 
total O abundance of about 8.45 dex, similar to the background 
average. 

It is interesting to compare the abundance distribution maps 
with the values that are derived from the extracted spectra of the 
proplyds. For the proplyds 170-337 and 170-334, chemical analysis 
reveals a complex abundance pattern because of the background 
contamination and the use of unreliable densities and temperatures. 
We have omitted the chemical abundances derived for these proplyds 
in Table |4| and, hence, focused on the results for 177-341, whose 
ionic abundances are presented in Table |3] We also show results 
considering the eflPects of opacity in this proplyd (see Taking 
into account the role of the collisional de-excitation mentioned in 
i3T]and 1421] it is clear that the flux of the [O ii] i3728 doublet in 
the observed spectra does not provide reliable abundances. This dou- 
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Table 4. Physical conditions for the background-subtracted spectra (c(I^)P = 0), the observed spectra (c(I^)P :» 1) and the adjacent background 
of the proplyds 170-337 and 170-334. 







170-337 






170-334 






Proplyd Emission 


Background 


Proplyd Emission 


Background 







»1 


Emission 





»1 


Emission 




0.76±0.22 


0.68±0.15 


0.72±0.15 


0.53±0.28 


0.70+0.15 


0.66+0.15 






Physical 


conditions 








«e(cm-3) 


1850±950 


3180±1600 


4080±2300 




6000: 


3760±2050 


[Fern] 




125900±39800 


5950±1280 






5380: 


diag. 


128800±29900 


53700±9950 




158500±36800 


24000±6100 




adopted^ 


128800±29900 


57900±9600 


5650±1170 


158500±36800 


24000±6100 


3760±2050 


re(K) [Nil] 


9470±1170 


8360±470 


10100±340 


10740±1660 


10100±770 


9590±580 


[Oin] 


9490+400 


8600+200 


8870+190 


10870+560 


10200+270 


8790+190 


diag. 


9430±490 


8560±240 




10770±690 


10060±240 





^ Determinations based on diagnostic curves (Fig.jTJ. 



^ Proplyd emission: weighted average of [Fein] and diag. Uq-Tq densities where possible. Background emission: weighted average of [Feiii] 
and [S ii] densities where possible. 



blet is not detected in the background- subtracted spectra of 177-341. 
Thus, the 0^/H+ and 0/H ratios have also been omitted in Table [3] 
We obtained an O^^ abundance of about 8.36 dex from CELs for 
the observed spectra, similar to the adjacent background value, 8.30 
dex, and higher than the abundance obtained from the 0^^/H+ ratio 
map (Fig. [S]?). Subtraction of the background emission yields an 
ratio of 8.42 dex for 177-341. The other ions show larger 
differences in abundance when comparing the observed, background 
and background- subtracted spectra of the proplyd. The largest dif- 
ferences with respect to the adjacent background are found for the 
and abundances of about 1 dex in the background- subtracted 
spectra and 0.5 dex in the observed one. These variations are related 
to the stronger dependence of the lines emitted by these ions on 
the assumed physical conditions. There is a further complication 
in the abundance calculations in proplyds because there may be 
a density stratification accross the objects. Rubin^ (^1989 ) pointed 
out that ionic abundances of heavy elements in the ionized gas de- 
rived from empirical techniques using CELs can be aff'ected by the 
presence of density changes along the line of sight. It is clear that 
this efl'ect would be even more severe when ionic abundances are 
calculated making use of emission lines with low critical densities. 
Hence, auroral lines should provide more reliable abundances than 
nebular lines due to their higher critical densities (see Tabl e [Tjl. In 
the particula r case of LV2, with a density of 7.9 x 10^ cm~^, Tsamis| 
|et al.| ( |201lt found that the O^^, and abundances derived from 
auroral lines were 10%, 30% and 50% higher than those from nebu- 
lar lines. Contrary to the results for LV2, we have found an excellent 
agreement between the N^, O^^ and abundances derived from 
nebular and auroral lines. For the N^/H+ and 0^^/H+ ratios, the 
diff'erences of using auroral and nebular lines are between 1 % and 
3% in both background- subtracted and observed spectra, which are 
equal to the diff'erences obtained in the adjacent background from 
both kinds of lines for these ions. For the abundances, we obtain 
an abundance difi'erence between auroral and nebular lines of about 
5% (0.02 dex) in the adjacent background. This diff'erence amounts 
to 35% (0.15 dex) in the background- subtracted spectra, but it is of 
the order of the uncertainty of the S^/H+ ratio. The observed spectra 
show an S^/H+ abundance from nebular lines that is 0.61 dex higher 
than that derived from auroral lines. This is because most of the 
[S ii] nebular line emission comes from the nebular background, 
while auroral lines are emitted from both proplyd and background 
components. It is encouraging that we obtain a very good agreement 
between nebular and auroral lines for the abundances of the ions 



above, which corroborates the validity of the background subtraction 
for 177-341 and the physical conditions determined for this proplyd. 



4.4 Chemical abundances from RLs 

For 177-341, we determined the He^ abundances using the Hei 
/14471, /15876, /16678 lines, which were weighted by 1:3:1 in the 
computation of the mean values according to their relative line fluxes 
(see Table [3]). We took the eflfective recombination coeflficients of 
Storey & Hummer] ( [1995| l for H i and those compute d by ,Porter| 
et al.|(|2005| t with the interpolation formulae provided by [Porter et al.| 
f 20071) for He i. The coUisional contribution was estimated from 
Sawey & Berrington| |l993 ) and Kingdon & Ferland ( 1995 ), and the 
optical depth in the triplet lines was derived from the computations 
of |Benjamin et al.| < [2002| ). The adopted physical conditions used to 
compute the abundances were from the diagnostic diagrams (see 
Table[3]l in both background- subtracted and observed spectra. For the 
background emission, we adopted an average value of low ionization 
zone, re([Nii]), and high ionization zone, re([Oin]), weighted by 
the relative ionic fractions 0^/0 and O^^/O, respectively. The results 
of Table [3] show a very similar He^/H+ ratio for the diff"erent dust 
opacity cases. 

The good signal-to-noise ratio of our spectral data allows us to 
detect and measure pure RLs of 11 and C n in most spaxels of the 
field and, therefore, also in the extracted spaxel of the proplyd 177- 
341. In Fig.|2]we show a section of the 177-341 spectrum around 
the lines of IMultiplet 1 of 11. RLs of heavy-elements present the 
advantage that their line flux ratio with respect to an H i line depends 
very weakly on the physical conditions, and thus their measurements, 
and consequently the abundances derived, are much less aff'ected by 
temperature fluctuations (see P eimbert|1967| l and, in our case, the 
presence of density inhomogeneities. 

The spectral resolution of our data allows to measure the blend 
of the brightest lines of IMultiplet 1 of 11: 11 /14649, /14651. Then, 
we used the empirical relations (1) and (4) obtained by 'Peimbert &J 
Peimbert ( 2005 ) to estimate the expected flux of all 11 Ml multiplet 
lines with respect to H5 and determine the O^^ abundances. Those 
expressions also correct for the NLTE eff'ects in the relative flux 
of each line of the multiplet, though these eflfects are expected to 
be rather small in the Orion Nebula because of its relatively high 
density, and probably negligible in denser proplyds. Combining the 
previous empirical relations, we can calculate the 0^^/H+ ratio from 
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Figure 9. Spatial distribution maps of C^^/H+ (a) and O^^/H^ (b) abundances derived from RLs and the map of ADF(0^^) (c). Abundances are presented in 
units of 12 + log(X+VH^). Ha contours are overplotted in all maps. RLs are not detected in the black areas at the upper right corner of the field (positions of 
proplyds 170-334 and 170-337). 



RLs as: 

^ ^ ^ QfefF(H/?) ^ 7(0 II Ml) 

H+ 4861 ^ aeff(Ml) ^ 1(11/3) ' ^ ^ 

where a (11/3) and Qfeff(Ml) are the effective recombination co- 
efficients for 11/3 and for the On Ml multiplet, respectively, and 
Ami = 4651.5 A is the representative mean wavelength of the whole 
multiplet. We can also determine the C^^/H^ ratio from the flux 
of C II /14267 RL using a similar relation. According to the two- 
zone scheme defined in ^4.3| re([Oiii]) was adopted to calculate 
the abundances of the two ions, 0^+ and C^^. We adopted the On 
effective recombination coefficients from |Storey| ( | 1 994| l, who carried 
out the calculations in the LS coupling scheme, and the C ii effective 
recombination coefficients from Davey et al. (2000 ). The typical un- 
certainties amount to 0.11 dex for both C^^/H^ and O^^/H^ ratios. 

The spatial distribution of the C^^ and 0^+ abundances from 
RLs are presented in Fig. [9] and the values obtained from the anal- 
ysis of the spaxel encompassing 177-341 are shown in Table|3] It 
should be noted that we have not detected RLs in the spaxels at and 
around the proplyds 170-337 and 170-334 (the black areas of Fig.[9]). 
In Fig.[9];a)the C^+ZH^ ratio map shows an almost featureless distri- 
bution reaching its highest value near the proplyd 170-334. However, 
this may be an artifact because at those spaxels the C ii /14267 line is 
affected by emission features of the neighboring proplyd. We found 
that the O^^/H^ ratio map obtained from RLs (Fig.[9p) has similar 
behaviour to that derived from CELs (Fig.[8]3), except for the spaxels 
(2 : 0, -3 : 0) that correspond to the position of the bow shock in front 
of 177-341, where the 0^+ abundance from RLs shows an enhance- 
ment that is not seen in the other maps (see ^4.5| >. This fact is not due 



to a low signal-to-noise ratio of the O ii RLs because they are rather 
bright at this zone. Finally, it is interesting to compare the results of 
the proplyd 177-341 on the maps with those derived from the ob- 
served spectra and included in Table[3] On the maps and for the same 
spaxel of 177-341 (see ^2.2\ , we obtain C^^ and 0^+ abundances of 
8.32+0. 1 1 and 8.42+0.1 1, respectively, which are very similar to the 
values derived from the observed spectra case of Table|3] 8.33 + 0. 1 1 
and 8.40 + 0.1 1, respectively. This is as expected because of the weak 
dependence on physical conditions of the abundances derived from 
RLs. Moreover, due to this weak dependence, the values obtained for 
the adjacent background emission (12 + \og(C^^ /H^) = 8.35 + 0.07 
and 12 + log(02+/H+) = 8.46 + 0.11) and those derived from the 
background-subtracted spectra (12 + log(C2+/H+) = 8.25 + 0.12 and 
12 + log(02+/H+) = 8.37 + 0.14) are very similar considering the 
relatively low signal-to-noise ratio of these RLs. Therefore, the dif- 
ferent abundance determinations from RLs, including also those of 
He^, remain rather constant and unaffected by the dramatic density 
contrast that we observe in Table [3] 

4.5 Abundance discrepancy factor of 0^+ 

By comparison of the O^^/H^ abundances determined from RLs 
and CELs we have computed and mapped the AD factor of O^^, 
ADF(0^^), adopting the following definition for this quantity: 

In Fig. [9|c), we can see that the ADF(0^^) map presents values 
between 0.10 and 0.35 dex. We do not compute the ADF(0^+) at 
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170-337 and 170-334 and their immediate surroundings because 
the O II RLs are not detected in those areas. Excluding the proplyd 
positions, we obtain an average ADF(0^^) of about 0.23 dex with 
a standard deviation of 0.06 dex for the whole map. This value 
is similar to previous determinations in other parts of the nebula 
( [Esteban et al.|2004[ [Mesa-Delgado et al.|2008| [20091 [20TT] ). As 
mentioned in the previous section, the spaxels at the position of 
the bowshock related to 177-341 show an increase in the flux of 
O II RLs but not in the [O iii] CELs. This produces an ADF(02+) of 
about 0.29 dex, while its immediate surroundings show values about 
0.22 dex. Deeper and more detailed spectroscopic observations are 
needed in order to confirm and explore this anomaly. 

From Table [3] we can see that the ADF(0^^) of the observed 
spectrum of 177-341 amounts to 0.04 + 0.11 dex, while in the 
ADF(0^^) map the value of that precise spaxel is about 0.30 dex. 
This dramatic diff'erence is related to the and bias aff'ecting the 
CELs when nebular contamination is not taken into account. The 
reasoning is the following: a) we use the densities of the /ie([S 11]) 
map, which are too low for the proplyds as discussed in g4.2.1[ for 
the calculations; b) the re([Oiii]) map derived using those densities 
produces artificially higher temperatures at the proplyds; and c) 
the O^^ abundances calculated from CELs in this way are lower 
than expected. In Table [3] we can also see that the ADF(0^^) of 
the background-subtracted spectra of 177-341 is -0.06 + 0.15 dex. 
Comparing the O^^ abundances determined for the difi'erent cases 
in Table [3] it can be seen that the change of the ADF(0^^) from 
one case to the other is mainly due to variations of the 0^^/H+ ratio 
derived from CELs. Those variations are larger than the quoted 
errors, in contrast to what happens for the abundances determined 
from RLs, which can be considered basically the same within the 
errors in the background- subtracted spectra and in the observed one. 
A similar result is found for the other ionic abundances obtained 
from RLs: He^ and C^^. Therefore, the variations of ratio 
derived from CELs may be produced by the temperature structure 
changes related to the presence of high-density ionized gas (and 
perhaps strong density stratification) in the proplyd. In contrast, 
abundances derived from RLs remain unaff'ected due to their weaker 
dependence on the physical conditions. Similar results were also 
found in the proplyd LV2 analysed by [Tsamis et al.[ ( [201 l| l and 
clearly show that the presence of high-density ionized gas is playing 
a role in the AD problem. 

Based on an IFU analysis of LV2 and the surrounding nebula, 
and taking into account theoretical considerations ( |Rubin| 1 989) [ Vie-[ 
[gas & Clegg[1994| t, [Tsamis et al.[p01l| l concluded that an important 
contribution to the AD problem in H 11 regions is the presence of 
dense regions (e.g., proplyds, globules or filaments). If these results 
can be generalized, the existence of a population of high-density, 
small and semi-ionized clumps/filaments, mixed with the nebular 
background gas of the Orion Nebula and randomly located along 
the line of sight, may be aff"ecting the classic [O iii] /14363//14959 
forbidden line temperature ratio within the observation aperture and, 
therefore, the gaseous abundances derived from CELs. In that case, 
CELs would yield lower limits to the mean Orion abundances, which 
would be more reliably determined from the faint metallic RLs. This 
last possibility is supported by Simon-Diaz & Stasiriska 1201 Ijl, who 
conclude that the stellar O abundances of the Orion OB Association 
are more consistent with gaseous ones determined from RLs (after 
correcting for the oxygen depletion onto dust grains). Further to be 
addressed is the question of how any such density inhomogeneities 
relate to the classic temperature fluctuation paradigm ( Peimbert 
[1967| l, which has been invoked as a possible solution to the AD 



problem. Follow-up studies on the link between the ADF(0^^) and 
temperature fluctuation parameter, f-, determinations are needed. 

Considering the results of LV2 ( (Tsamis et al.|201l] l and 177- 
341 (this work), we can conclude that these unknown clumps, similar 
in properties to the proplyds, need not be strong metallic RL emit- 
ters relative to the background emission of the Orion Nebula. This 
is supported by the fact that both studied proplyds show line flux 
ratios with respect to H/? of the RLs He i, 11 and C 11 that are sim- 
ilar to those of the nebular background (see Table [3] for 177-341). 
This line of reasoning is difl'erent from that of previous scenarios 
in which density inhomogeneities in H 11 regions were assumed to 
be hydrogen-deficient with respect to the host nebula and strong 
RL emitters on account of their low electron temperatures ( [Tsamis[ 
[& Pequignot[[2005[ [Stasinska et al.|2007| l. Another important is- 
sue is the size of the posited clumps and their spatial distribution, 
which would determine the total amount of ionizing radiation that 
they receive and the spectrum that they emit. In this sense, previous 
studies performed in the Orion Nebula (e.g. [Rubin et al.[2003][Mesa-| 
"Delgado et al.'2008) have not found any evidence for clumps, other 
than the known proplyds. The work carried out by [Rubin et al.[p003] t 
based on long-slit spectroscopy obtained with the HST attained a 
spatial resolution of 0'.'5, while ,Mesa-Delgado et al.[ ( p008| ) studied 
the spatial distribution of physical conditions and abundances at V!2 
from ground-based observations. Neither study detected significant 
small-scale density variations in the nebular background. On the 
one hand, in [Mesa-Delgado et al. ( 2008 ) the measurements were 
made of We([S n]), which can be affected by collisional de-excitation 
in the presence of high-density gas. On the other hand, [Rubin et al.[ 
( [2003] ) used a low density limit for the temperature determinations. 
So, it is clear that deeper spectroscopic studies using the HST at 
the highest spatial resolution (e.g., with WFC3 in the UV/far-UV), 
along with the use of appropriate density indicators are needed to 
explore the astrophysics of small scale "dumpiness" in the Orion 
Nebula. Finally, it should be noted that if these kinds of clumps are 
real, or simply there are important density variations in the nebular 
background, they would severely modify the CEL diagnostics. The 
existence of dense clumps within more distant Galactic and extra- 
galactic Hii regions can complicate the interpretation of emission 
line diagnostics with important ramifications for chemical evolution 
studies. 



5 DUST EFFECTS ON NEBULAR PROPERTIES OF THE 
PROPLYD 177-341 

Dust particles in the inner dense neutral zones of the proplyds can 
absorb and scatter the emission from the surrounding nebula, which 
complicates the background subtraction process and the determina- 
tion of the intrinsic proplyd emission. However, the dust properties 
of the proplyds are not well known. The only existing study is that 
of [Rost et al.[ ( [2008] l, who attempted to model near-infrared polar- 
ization mapping of 177-341, but found that simple models do not 
reproduce the observations. We therefore adopt a simplified empiri- 
cal procedure, similar to that presented in Figure 3a of jHenney &| 
0'Dell| ( jl999| l. The total dereddened flux for a given wavelength, /], 
measured at a proplyd position can be expressed schematically as: 

= /J, (3) 

where represents the nebular background emission, which is 
attenuated by an effective opacity proplyd optical depth, r^, and /J 
is the intrinsic proplyd emission. 

The effective optical depth will be significantly smaller than 
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the true absorption optical depth of the proplyd due to three effects. 
First, some portion of the nebular emission may arise from in front 
of the proplyd and so will not be absorbed. Second, the neutral core 
of the proplyd is small compared with the angular resolution of the 
observations, which produces "beam dilution" of the absorption. 
Third, scattering of nebular light by dust in the proplyd will partially 
"fill in" the absorption. This last eff'ect will depend on the geometrical 
relation between the proplyd and the bulk of the emitting gas in the 
nebula (as well as on the dust properties). For 177-341, the scattered 
contribution might be relatively high because we clearly see the 
proplyd projected against a particularly faint part of the nebula, 
whereas there are much brighter parts oflF to the side. In principle, 
we could constrain this by analysing the scattered continuum, but 
it requires further investigation, which is beyond the scope of this 
article. 

In the absence of clear observational constraints, we adopt 
a standard reddening law, f(A,Ry), parametrized by the total-to- 



Table 5. Input parameters for example physical model of 177-341 



selective extinction, Ry (Cardelli et al. 1989|, which will be a func 



tion of the grain optical properties and the relative importance of 
scattering. The proplyd opacity at any wavelength A is then related 
to the H/? extinction coefficient by 

p ciH/sy 



log 6 



■ (i+/(^,0) 



and the absolute dereddened flux of a proplyd is given by 

jP ^ Jt _ ^b^Q-c(H^)P(l+/y,J!P)) 



(4) 



(5) 



Assuming that the foregroud extinction, c(H/?)^°% is the same for the 
proplyd and adjacent nebular emission, then the raw observed fluxes 
(F^ = he'^^ ) follow a parallel relation: 



(6) 



This assumption is justified for the proplyd 177-341 since it can 
be seen from Table |3] that the c(H/?)^°^ coefficients derived from 
the background- subtracted, observed, and background spectra are 
identical within their respective uncertainties. 

It is also interesting to note that the cases presented in ^2.2| 
-and analysed along this paper for the diflPerent proplyds- can be 
easily obtained from Eq. 6. We can see that the transparent case, 
c(H/?)P = 0, is just a simple subtraction of the nebular spectra from 
the total observed spectra, where both nebular and proplyd emissions 
are mixed. On the contrary, in the opaque case, the proplyd spectra 
directly corresponds to the observed one. An interesting point of 
these two cases is that neither of them depend on the dust properties 
of the proplyd. 

These relations allow us to go beyond the approximations used 
in the previous sections, in which the proplyd was assumed to be 
either fully opaque (c(H/?)p » 1) or fully transparent (c(H/5)p = 0), 
and to explore the eflFects of partial transparency on the derived 
physical conditions and chemical abundances of 177-341. The final 
results are presented in Table |3] As well as for the transparent and 
opaque cases, we show diagnostics derived assuming c(H/?)p = 0.05 
and 0. 1 dex, which are similar to the values determined by Henney &, 
|0'Dell| ( [T999l l for this proplyd (0.04 to 0.08 dex). This pair of values 
is combined with three diff'erent values of Ry - 3.1, 5.5 and 7, using 
the f(A) reddening law of Blagrave et al.| (|2007|). Once we have 
estimated the proplyd fluxes, F^, for the diff'erent cases, we proceed 
with the nebular analysis calculations: c(^5)^°^ physical conditions 
and chemical abundances. As we see in Table [3] the standard ISM 
value for 7?^ of 3.1, produces c(H/?)^°^ coeflficients that are similar 
between the proplyd and background emissions, while higher Ry 
would imply lower c(H^?)^°^for the proplyd. This indicates that the 



Stellar spectrum: 



I Simon-Diaz et al, 



2006 
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Ionization front radius: 
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= 2.04 X 10^ Lq 

Oh 
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He 10.98, C 8.41, N 7.85, 8.30, 
Ne 7.56, S 6.98, Ar 6.26, Fe 5.65 
Standard Orion {Baldwin et al. 199l) 



dust properties in the proplyd 177-341 might be better represented 
by Ry = 3.1, though the uncertainties are too large to reach a solid 
conclusion. 

With respect to the physical conditions, we can note that the 
determinations are the most aflfected by the proplyd extinction. For 
R^ = 3.1, the adopted density falls from 3.8 x 10^ to 2.4 x 10^ cm'^ 
when c(H/?)P increases from to 0.1 dex. A similar contrast is found 
for the other two values of Ry. However, the electron temperatures 
do not show large diflFerences. Indeed, they slightly decrease for 
higher/?^ ratios. For the same c(H/?)p cases and Ry = 3.1, re([Nn]) 
and re([Oin]) show variations of 200 and 130 K, respectively. The 
derived chemical abundances are in turn aflfected by the changes 
in the physical conditions. The ionic abundances of N^, Ne^^, S^ 
and S^^ are the most aflPected by c(H/?)p variations because their 
abundance determinations have a stronger density dependence than 
the other ions. The largest diflFerences are found for S^/H+ and 
N^/H^ ratios: 0.24 and 0.19 dex, respectively when c(H/5)p varies 
from to 0.1, which is larger than the quoted uncertainties. Due to 
the much weaker dependence on the assumed physical conditions, 
ionic abundances of He^, C^^ and O^^ determined from RLs remain 
unaflfected. 



6 A PHYSICAL MODEL OF 177-341 

As an alternative to the purely empirical analysis presented in the 
previous sections, a diflFerent approach to analysing the emission 
spectrum of the proplyds is through the construction of physical 
models that combine a priori simulations of radiative transfer, hy- 
drodynamics, and atomic physics in order to predict the density, 
temperature, and ionization structure of the the photoevaporation 
flow through the proplyd ionization front. Such models have previ- 
ously been applied to the ensemble properties of large numbers of 
proplyds ( [Johnstone, HoUenbach & Bally|1998[ [Henney & Arthur] 
1998|l and in de tail to individual objects such as 177-341 fHenneyl 
&0'Dell|1999t , LV2 ( [Henney et al. r2002 ), and LVl (Graham et al] 
2002| ). We have improved on these models in two significant ways. 
First, whereas published models have considered only emission 
from regions where hydrogen is fully ionized and the flow is super- 
sonic, we now use a detailed analytic model of gas acceleration in 
the ionization front ( [Henney et al.|2005} to extend the treatment to 
cover partially ionized emission zones where the gas moves sub- 
sonically. Second, whereas published models used ad hoc fitting 
functions to the emissivity structure, specifically tailored to only the 
brightest emission lines, we now use the plasma microphysics code 
Cloudy ( Ferland et al. 1998 ) to self-consistently calculate the full 
physical structure and emission spectrum of the proplyd flow. 

The input parameters for the model (values given in Table|5]) 
are the radius of the ionization front at the proplyd cusp (assumed 
hemispherical), the intensity and spectral shape of the illuminating 
stellar radiation, and the composition (gas-phase elemental mix and 
dust grain populations) of the proplyd material. For most of these 
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Figure 10. Emission structure in the Hq-Tq plane for the example physical model of 177-341. (a) Positive-colored image of the three emission lines: [S ii] /16731 
(red), [Nil] /16583 (green), and [Oni] /15007 (blue), with brightness proportional to the fraction of the model luminosity in each line, (b) Negative-colored image 
of the same data overlaid with the observational diagnostic curves from Fig.|7]([S ii] in red, [Nii] in green, [Om] in blue). 
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Figure 11. Comparison of the model and observed line ratios for the proplyd 177-341. The model parameters are given in Table[5]and the emission structure of 
the model is presented in Fig.FTol Observed intensities correspond to the background-subtracted values (c(I^)P = 0) given in column 4 of Table 2. 



parameters, we have taken values from the literature, whereas the 
gas-phase abundances and incident ionizing flux have been adjusted 
slightly in an attempt to reproduce the observed emission line in- 
tensities. The stellar spectrum is that determined for the ionizing 
star 6^ Ori C by spectroscopic analysis ( |Sim6n-Dfaz et al. |2006| l. 
The ionization front radius is the value determined by fitting to HST 
emission line images ( Henney & O'Dell 1999 ), while the adopted 
ionizing flux at the proplyd position corresponds to a physical sepa- 



ration of the proplyd from the ionizing star of2.13xl0^^ cm if there 
is no intervening absorption. Given the observed angular separation 
of 25784 (Bally et al. 1998 |, and assumin g a distance to the Orion 
Nebula of 440 pc ( O' Dell & Henney|2008] t, the projected separation 
is l.VOx 10^^ cm, implying an inclination angle of ^ 55°. The diff'use 
radiation field and the proplyd tail are ignored in this model, since 
the observational aperture (see ^2.2| l only covers the head of the 
proplyd. 
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The basic structure of the model is of a slow ionized wind that 
expands spherically away from the ionization front at the bright 
proplyd cusp (see figure 1 of |Henney & Ar thur 1998 ). The outer 
layers of this externally ionized wind have the highest ionization 
parameter because they are relatively low density (< 1000 cm~^) 
and see the full ionizing flux from the Trapezium stars. The inner 
layers of the wind are denser and see a reduced flux because of 
absorption in the overlying layers. This gives a lower ionization 
parameter and also a higher temperature due to radiation hardening 
and collisional quenching of coolant lines. The temperature reaches 
a maximum value of ~ 11000 K just outside the ionization front, 
while the electron density peaks inside the front itself. For the in- 
nermost layers, hydrogen becomes increasingly neutral and both 
the temperature and electron density drop rapidly. Thus, each radius 
within the proplyd has its own combination of density, temperature, 
and ionization state. It is the correlations between these quantities 
(as illustrated in Fig. [T0| that determine the optical emission line 
spectrum of the proplyd (Fig.pT|). 

Fig.[ToJa) shows the emission structure of [Sii] /16731 (red), 
[Nil] A65S3 (green), and [Oiii] /15007 (blue). The low-density outer 
zones of the model are highly ionized, emitting principally in [O iii], 
as seen in the blue horizontal branch at ~ 9700 K in the upper 
left part of the figure. As one approaches the proplyd ionization 
front (direction of increasing density), the [N ii] and then [S ii] emis- 
sions become relatively stronger (color changes to yellow) and the 
temperature increases, reaching a maximum of ~ 1 1000 K. In the 
ionization front itself, where the gas is partially neutral, the [S ii] 
emission comes to dominate over [N ii], giving rise to the orange-red 
branch that curves down and to the left. Fig.fTofb) shows a negative 



Tola) overlaid with the 
a) ([Sii] in red, [Nii] in 



representation of the same data as in Fig. 
observational diagnostic curves from Fig.|7 j 
green, [O iii] in blue). The white star shows the solution for («e, T^) 
obtained in ^4. 2 [ under the assumption that the emission in all three 
ions is co-extensive at a single density and temperature. This solu- 
tion lies well outside the locus of densities and temperatures seen in 
the model, even though the model does reproduce the observed line 
ratios (see Fig.[TT]». Instead, the model satisfies the constraints of the 
three diagnostic curves separately in three diff'erent regions, which 
are indicated schematically by the colored stars. Note, however that 
even within the emission region of each ion there is considerable 
variation in the physical conditions. Independent evidence that the 
higher proplyd temperatures deduced from the physical model may 
be correct comes from the temperature map of Garci a-Diaz et al.| 
( [2OO8I 1, who found temperature peaks of 10,000 to 12,000 K at the 
positions of several proplyds, including 177-341 (see their figure 
21). Those authors used a linewidth method instead of the tradi- 
tional line ratio methods, and so their results should be insensitive 
to uncertainties in the electron density. 

Fig. [TT] shows that the model generally reproduces very well 
the observed relative line intensities for the intrinsic spectra of 177- 
341 with the notable exception of the [O i] /16300 line, which has a 
predicted flux of only half the observed value. This discrepancy may 
be explained if a significant fraction of the observed [O i] emission 
comes from dissociation of OH at the surface of the proplyd' s cir- 
cumstellar disc fStorzer & Hollenbach|1998| l, which is not included 
in our model. Another factor that may increase such discrepancy is 
that the [O i] emission in the observed and background spectra has 
not been corrected from telluric emission contamination, though the 
subtraction process should minimize it in the intrinsic spectra of the 
proplyd. The only other disagreement is with the [S iii] /16312 line, 
which is roughly 25% too high in the model (roughly 1.5 times the 
observational uncertainty). This may be an indication that the details 



of the stellar spectrum around 23 eV are not correctly modeled, 
although the disagreement is only marginal. 

Given that the stellar spectrum is very well constrained, the 
only way to vary the model temperature by a significant amount is by 
changing the gas-phase metal abundances and thus the gas cooling 
rate. In particular, the temperature in the highly ionized outer regions 
of the proplyd is very sensitive to the oxygen abundance since the 
[O III] /15007 line is the major coolant there. In order to reproduce 
the observed line intensities (particularly the [O iii] /14363//15007 
ratio), we find it necessary to reduce the oxygen abundance to less 
than 50% of the value that has been measured for the surrounding 
nebula (e.g., Esteban et al.||200 4), which leads to a temperature 
of around 9500 K in the [Oiii] emission zone. Using the standard 
Orion gas-phase abundances gives temperatures of about 8500 K 
for the same zone and predicts [O iii]/B/3 ratios that are considerably 
higher than observed for the nebular [O iii] /14959, A5001 lines, but 
lower than the observed for [O iii] /14363 auroral line. Using the even 
higher abundances derived by |Tsamis et al.| ( [2011| l for the proplyd 
LV2 implies temperatures as low as 7700 K and exacerbates the 
disagreement with the observations still further. For comparison, 
both varying the stellar eff"ective temperature by +1000 K (the un- 
certainty in the analysis of [Simon-Diaz et al.|2006| l or varying the 
ionizing flux by a factor of two (corresponding to the uncertainty in 
the inclination angle of the proplyd) only change the gas temperature 
by < 200 K, while varying the stellar atmosphere model between 
WMBasic ( Pauldrach et al.|2001| l and TLUSTY ( |Lanz & Hubeny] 
i2003j ) produces a slightly larger change of ^ 400 K. 

Despite the promising success of the model in reproducing the 
observed line ratios of the intrinsic spectra of 177-341, further work 
is necessary before a definitive statement can be made about the gas- 
phase abundances in the proplyd. The lack of observations of [0 11] 
lines, together with the uncertain nature of the [O i] emission (see 
above) means that the oxygen abundance hinges on observations of 
a single ion stage. It is vital to test the model against other proplyds, 
particularly those such as LV2 ( Tsami s et al. 2011 , Tsamis & Walsh] 
|20TT] l where the [On] /13726, ^3729 doublet is observed. At the 
same time, the models need to be simultaneously constrained not 
just by the line ratios, but also by sub-arcsecond HST imaging 
(Bally et al. 1998, O'Dell 199 8) and velocit y profiles from high- 
resolution spectroscopy (Henney & O ' Dell]] 1 999 [ [Shuping et al.| 
|2003| ). It is also necessary to examine in greater detail the potential 
influence of internal extinction within the proplyd. This depends on 
the properties of the dust that is entrained in the photoevaporation 
flow, which can be constrained by comparison with observations 
of thermal mid-infrared emission ( [Garcfa-Arredondo et al.|200H 



[Smith et al.|2005l|Robberto et al.|2005]|Shuping et al.|2006| ). 



7 CONCLUSIONS 

In this paper we used integral field spectroscopy in the spectral range 
3500-7200 A to explore the eff'ect of the presence of protoplanetary 
discs (proplyds) on the controversial abundance discrepancy (AD) 
problem. Our field is located to the south of the Trapezium Cluster 
and contains difi'erent morphological structures. Among them, the 
most prominent structures are the proplyds 177-341, 170-337 and 
170-334 and the bowshock in front of 177-341 and the high- velocity 
jet HH 514 emerging from 170-337. 

We analyse spatial distribution maps of several nebular proper- 
ties: emission line fluxes, extinction coefi&cient, electron densities 
and temperatures, chemical abundances from collisionally excited 
lines (CELs) and recombination lines (RLs), and the AD factor of 
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O^^, ADF(0^^). We found that collisional de-excitation due to high 
densities is the main mechanism to explain the behaviour observed 
in the spatial distributions of the nebular properties at the position 
of the proplyds. Analysing the collisional de-excitation effects, we 
estimate that proplyds 177-341 and 170-334 should have a density 
between 9 x 10^ and 7x10^ cm~^, the critical densities of the depar- 
ture levels of the [Nn] and [Om] nebular lines, respectively. For the 
proplyd 170-334, we estimate that it should present a density higher 
than 7 X 10^ cm~^, the critical density of [Om] nebular lines, based 
on the intense enhancement observed for this proplyd in the spatial 
distribution of the [O ni] auroral line. 

In this study we also obtain the intrinsic spectra of the pro- 
plyds 177-341, 170-337 and 170-334 by subtracting the nebular 
background emission. From a detailed analysis, we derive densities, 
temperatures and chemical abundances using empirical diagnostics 
for the three proplyds. Due to their smaller sizes, the diagnostics re- 
veals that proplyds 170-337 and 170-334 present some background 
contamination, so the physical conditions and chemical abundances 
derived for them may not be correct. 

We detect O ii and C ii RLs in the spectrum of the proplyd 177- 
341. The density determined in its background- subtracted spectrum 
is about 3.8 X 10^ cm~^, which is much higher than the density of the 
nebular background of the Orion Nebula, 3800 cm~^. We determine 
electron temperatures from [N ii] and [O m] line ratios finding that 
the striking spikes observed in the spatial temperature distribution 
maps disappear once the correct density value for the proplyds is 
used.. We obtain the chemical abundances of an extensive set of ions 
from CELs: N+, 0^+, Ne^+, S+, S^+, Ar^^ and Fe^^. We also deter- 
mine ionic abundances from RLs for He^, O^^ and C^^ and observe 
that ionic abundances from CELs show important difi'erences when 
comparing the observed spectrum and the intrinsic one of 177-341, 
while those derived from RLs remains similar. We found that ionic 
abundances of N^, O^^ and derived from nebular and auroral 
lines are very consistent with one another when using the correct 
high densities. 

We investigate the eff'ects on our results of internal extinction 
by dust in the proplyd's neutral core. We find that the resultant uncer- 
tainty in the background subtraction process may have a significant 
impact on the derived density and ionic abundances (such as S^/H+ 
and N^/H+) that are sensitive to this. 

As an alternative to the empirical analysis we construct a physi- 
cal model for the proplyd 177-341, obtaining the ionization structure 
of the proplyd as well as its temperature and density. We find an 
excellent agreement between the observed line ratios of the intrinsic 
spectra and the predicted ones from the model, even though the 
predicted physical conditions for the well-resolved S^, and O^^ 
zones diff'er from those obtained in the one-zone scheme applied in 
the empirical analysis. Regarding the final gas-phase abundance in 
the proplyd, further work is needed to definitively determine them. 

Finally, comparing the 0^^/H+ ratio derived from CELs and 
RLs in the proplyd 177-341, we find that the ADF(02+) tends to 
in the observed spectra and in the intrinsic spectra of this proplyd. 
A similar result was previously found in the analysis of the proplyd 
LV2 performed by |Tsamis et al.| ( [20TT| ). Though our determinations 
show relatively high uncertainties, we conclude that high-density 
gas (in the form of proplyds, globules or unseen clumps/filaments 
acting at small spatial scales) is playing a major role in the AD 
problem. Based on the results on 177-341 and LV2, we discuss 
the AD problem in the context of density inhomogeneities finding 
that the proposal of |Viegas & Cleggl ( |1994) seems to explain the 
behaviour observed in the proplyds. In a scenario where small, high- 
density and semi-ionized clumps/filaments are mixed with difi'use 



gas in the observation aperture, the classical CEL abundance anal- 
ysis can be afi'ected if collisional de-excitation of certain emission 
line diagnostics is not accounted for. RLs would not be aff'ected and, 
therefore, they should reliably yield the chemical abundances in the 
target field. Judging by the results of the proplyd IFU analyses thus 
far, these clumps need not be strongly hydrogen-deficient, unlike 
those posited in previous scenarios to explain the AD problem in 
H II regions ( Tsamis & Pequignot 2005 , Stasihska et al. 2007 ) and 
planetary nebulae ( (Harrington & Feibelman|1984^|Liu et al.|2000[ 
ITsamis et al.]2008] ). 
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